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Statement  of  the  Problem  Studied 

A  new  in  situ  method  for  characterization  of  unsaturated  soil  hydraulic  properties  was 
studied.  Knowledge  of  these  soil  properties  is  necessary  for  predicting  the  movement  of  water 
and  water-borne  contaminants  from  surface  and  near-surface  sources  to  groundwater,  design  of 
bioremediation  systems,  prediction  of  infiltration  capacity,  and  determining  the  susceptibility  for 
landslides  and  floods.  Many  existing  approaches  for  measurement  of  hydraulic  properties  in 
unsaturated  soils  are  limited  to  near-surface  application  and  require  augering  a  borehole,  which 
may  result  in  contaminated  soil  being  brought  to  the  surface.  In  this  project,  a  new  cone 
penetrometer  test  (CPT)  method  for  determining  the  hydraulic  properties  of  unsaturated  soils  was 
developed  and  tested.  While  CPT  methods  exist  for  measurement  of  the  hydraulic  conductivity 
of  saturated  soil,  no  methods  were  previously  available  for  determination  of  the  soil-moisture 
characteristic  curve  and  hydraulic  conductivity  function  of  imsaturated  soils.  This  work  is 
important  because  CPT  methods  for  site  exploration  are  preferred  at  contaminated  sites.  CPT 
methods  do  not  produce  soil  cuttings  at  the  surface,  and  can  provide  more  rapid  characterization 
of  the  subsurface  than  traditional  methods.  Given  the  potential  of  this  new  CPT  method  for 
expedited  site  characterization,  this  work  directly  supports  the  Army’s  needs  improved  methods 
for  characterizing  and  remediating  its  contaminated  sites. 


Research  Objectives 

The  overall  goal  of  this  project  was  to  develop  a  CPT-based  methodology  for 
investigating  the  fundamental  relationships  between  measurable  transient  unsaturated  flow 
responses  and  soil  hydraulic  parameters.  Specific  objectives  of  this  work  included  1)  designing  a 
modified  cone  penetrometer  (called  a  cone  permeameter)  to  inject  water  into  to  subsurface  and 
measure  pressure  head  responses,  2)  creating  an  inverse  methodology  for  estimating  hydraulic 
properties  of  unsaturated  soils  from  cone  permeameter  data,  3)  demonstrating  the  performance  of 
a  prototype  via  carefully  controlled  full-scale  laboratory  tests,  and  4)  testing  the  methodology 
and  instrumentation  in  the  field. 
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Background 


The  cone  permeameter  was  designed  to  determine  the  soil  water  retention  curve,  d(h),  the 
relationship  between  the  water  content  and  the  pressure  head  of  an  unsaturated  soil,  and  the 
unsaturated  hydraulic  conductivity  function  of  the  soil,  expressed  in  terms  of  pressure  head  or 
water  content,  K(h)  or  K(9).  Knowledge  of  these  curves  is  required  for  accurate  numerical 
modeling  of  variably  saturated  flow.  The  van  Genuchten  (1980)  analytical  functions  are 
commonly  used  to  describe  these  hydraulic  properties: 

1 

e,-e,  (i+|a/j|”)"' 

e,  =  i,  h>Q 

(1) 

K{h)  =  K.e'  [i  -  (i  -  ey ^ f  <  0 

K{h)  =  K„  h>{i 

(2) 

where  dg  is  the  effective  water  content  (L^L'^),  Or  and  0^  are  the  residual  and  saturated  water 
contents  (L^L‘^)  respectively,  a  (L‘'),  n  and  m  (=l-l/«)  are  empirical  fitting  parameters  (-),  Kg  is 
the  saturated  hydraulic  conductivity  (LT')  and  /  is  the  pore  cormectivity  parameter  (-).  The  pore 
connectivity  parameter,  /,  is  commonly  set  equal  to  0.5  (Mualem,  1976). 

To  develop  a  CPT  tool  and  methodology  for  determining  the  soil  hydraulic  properties 
K(h)  and  6(h)  for  vadose  zone  soils,  the  cone  permeameter  was  designed  to  inject  water  into  the 
soil  under  constant  pressure  through  a  screen  and  to  measure  movement  of  the  wetting  front  with 
tensiometers  above  the  screen  (Figure  1)  (Leonard,  1997;  Gribb  et  al.,  1998).  The  parameters 
describing  the  K(h)  and  d(h)  relationships  are  determined  via  analysis  of  cumulative  inflow 
volume  and  pressure  head  readings  using  an  inverse  solution  technique.  To  accomplish  the 
inversion,  a  cone  permeameter  test  is  numerically  simulated  with  the  radially  symmetric  form  of 
Richards  equation  and  the  van  Genuchten  (1980)  functions  of  K(h)  and  d(h): 


dt 


dh 

d 

f  dh  (W 

rK  — 

+  — 

K 

—  +  J 

dr  _ 

dz 

(dz  )_ 

(3) 


Then  an  objective  function,  0,  expressing  the  differences  between  flow  responses  measured  with 
the  cone  permeameter  and  those  predicted  using  the  numerical  model  with  parameterized  soil 
hydraulic  properties,  is  minimized: 

^iP,q,p)  =  ^vj  '^w,j[qj{x,t)-qj{x,t„b  )p  Pj{e„b  )Y  (4) 

j-\  i=i  7-1  /-I 
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where  the  first  term  on  the  right  side  represents  deviations  between  measured  and  predicted 
space-time  variables  (e.g.,  pressure  heads  or  moisture  contents  at  different  locations  and/or  times, 
or  the  cumulative  infiltration  rate  versus  time).  In  this  term,  mq  is  the  number  of  different  sets  of 
measurements,  and  riqj  is  the  number  of  measurements  in  a  particular  measurement  set.  Specific 
measurements  at  time  tj  for  the  yth  measurement  set  at  location  x(r,  z)  are  represented  by  qj*{x, 
tj),  gj(x,  ti,  b)  are  the  corresponding  model  predictions  for  the  vector  of  optimized  parameters  b 
(e-g-5  Gv’  Kg),  and  vj  and  w/y  are  weights  associated  with  a  particular  measurement 

set  or  point,  respectively.  The  weighting  factor,  vy,  is  given  by  the  inverse  of  the  number  of 
measurements  multiplied  by  the  variance  of  those  observations,  and  wjj  is  equal  to  1  in  this 
work.  The  second  term  on  the  right  represents  differences  between  independently  measured 
pj*(6i)  and  predicted  pj(6i,  b)  soil  hydraulic  properties  (e.g.,  d(h),  K{ff)  or  K(h)  data),  while  the 
terms  mp,  ripj,  vj  and  w/y  have  similar  meanings  as  for  the  first  term,  but  are  now  for  the  soil 
hydraulic  properties.  Minimization  of  the  objective  function  0  is  accomplished  by  using  the 
Levenberg-Marquardt  nonlinear  minimization  method  (Marquardt,  1963),  which  combines  the 
Newton  and  steepest  descent  methods. 


Figure  1 .  Prototype  cone  permeameter. 
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Summary  of  the  Most  Important  Results 

The  project  objectives  described  above  were  met  via  a  fully  integrated  numerical  and 
experimental  research  program.  Numerical  studies  were  undertaken  to  investigate  the  effects  of 
varying  hydraulic  soil  properties  and  test  conditions  on  observed  flow  responses.  A  fully 
functional  prototype  cone  permeameter  and  data  acquisition  system  were  designed  and 
fabricated.  Full-scale  soil  tank  experiments  were  performed  to  assess  the  validity  of  the 
numerical  model  predictions  and  the  applicability  and  ease  of  use  of  the  prototype  in  sandy  soil 
(Gribb  et  al.,  1998;  Kodesova  et  al.,  1998a, b,  1999a;  Simunek  et  al.,  1999)  (Figure  2).  Finally,  a 
field  demonstration  of  the  prototype  was  completed  (Kodesova  et  al.,  1999b).  The  information 
gained  on  unsaturated  flow  processes  through  this  work  provides  the  basis  for  future  work  with 
other  soil  types,  and  investigations  into  the  effects  of  soil  anisotropy  and  disturbance  due  to  cone 
penetration.  Specific  results  include  are  as  follows: 

•  We  modified  a  standard  size  cone  penetrometer  housing  to  allow  injection  of  water  into  the 
subsurface  through  a  screened  section  at  a  known  pressure  and  measure  the  pore  water 
pressure  at  two  locations  with  tensiometer  assemblies  (shown  in  Figure  1).  We  developed 
inverse  solution  methods  to  predict  the  soil  hydraulic  properties  from  data  obtained  during 
transient  flow  tests  and  objectives  functions  of  the  form  shown  in  Eq.  4  (Gribb  et  al,  1998; 
Kodesova  et  al,  1998b,  1999a). 


Figure  2.  Cone  permeameter  set-up  for  laboratory  aquifer  and  field  testing. 
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•  We  tested  the  cone  permeameter  in  a  laboratory  aquifer  and  compared  the  results  of  inversion 
to  results  from  standard  methods  including  in  situ  slug  tests,  and  laboratory  constant  head 
permeability,  pressure  plate,  capillary  rise  and  hanging  column  tests.  The  cone  permeameter 
results  compared  favorably  with  results  from  the  standard  methods,  suggesting  that  this  tool 
can  be  used  to  estimate  K(h)  and  6(h)  in  sandy  soil  (Gribb  et  al.,  1998;  Kodesova  et  al., 
1998a,b,  1999a). 

•  We  then  tested  the  cone  permeameter  at  three  sites  in  a  sandy,  relatively  homogeneous  soil 
(Figure  3).  We  compared  the  inversion  results  to  results  from  standard  methods  including 
falling  head,  pressure  plate,  capillary  rise,  and  multi-step  outflow  laboratory  tests  and  in  situ 
Guelph  permeameter  tests.  The  cone  permeameter  results  compared  favorably  with  results 
from  the  standard  methods,  suggesting  that  this  tool  can  be  used  to  estimate  K(h)  and  6(h)  in 
sandy  in-situ  soils.  Figure  4  shows  test  data  (cumulative  inflow  voltime  and  pressure  head 
readings)  and  simulated  results  for  a  single  test.  The  estimated  6(h)  curves  \vith  respect  to 
those  obtained  using  the  standard  methods  are  shown  in  Figure  5,  and  the  saturated  hydraulic 
conductivity  values  obtained  from  inversion  of  test  data  from  one  location  are  shown  in  Table 
1 .  Excellent  agreement  between  the  hydraulic  conductivity  values  was  obtained  (Kodesova  et 
al.,  1999b). 

•  We  analyzed  the  usability/practicality  of  the  prototype  and  test  method  and  modified  the 
procedure  and  equipment  accordingly.  Analysis  of  the  procedure  included  the  impact  of  cone 
permeameter  placement  (pushed  vs.  buried)  and  the  impact  of  the  initial  soil  moisture  content 
on  estimates  of  the  soil  hydraulic  properties.  We  found  that  the  effect  of  soil  densification 
caused  by  direct  push  was  minimal,  but  resulted  in  slightly  lower  values  of  6s  and  Ks  than 
when  the  prototype  was  buried  prior  to  testing  (Kodesova,  et  al.,  1998a,b)  and  that  additional 
inputs  of  initial  soil  moisture  content  are  helpful  for  more  precise  determination  of  the 
hydraulic  parameters  usually  investigated  {6s,  6r,  a,  n,  and  Ks)  (Kodesova,  et  al.,  1998b, 
1999a).  The  equipment  modifications  included  modifying  the  screen  to  prevent  clogging, 
improving  the  method  of  assembly  for  the  tensiometer  sections,  and  devising  a  method  to 
insert  the  permeameter  into  the  soil  to  minimize  disturbance. 

•  We  also  examined  the  impact  of  varying  the  number  of  hydraulic  parameters  to  be  estimated 
and  found  that  optimization  of  additional  parameters  (/,  the  pore  connectivity  factor,  and  k^, 
an  anisotropy  factor)  improved  the  fit  of  measured  data  (Simunek  et  al.,  1999;  Kode§ovd,  et 
al.,  1999a).  We  also  showed  that  we  can  find  the  wetting  branch  of  soil  hydraulic 
characteristics  from  the  main  wetting  part  of  experiment,  as  well  as  simultaneously  evaluate 
the  wetting  and  drying  curves  via  analysis  of  both  wetting  and  redistribution  parts  of  a  cone 
permeameter  test  (Simunek  et  al.,  1999;  Kodesova,  et  al.,  1999a). 

•  In  addition  to  using  this  method  to  predict  soil  hydraulic  properties  of  unsaturated  soils,  we 
examined  the  feasibility  of  using  the  cone  permeameter  and  a  simplified  analysis  method  to 
predict  the  saturated  hydraulic  conductivity  rapidly,  while  still  in  the  field  (currently,  the 
inverse  solution  is  not  suited  to  analysis  in  real  time).  We  found  that  the  log-transformation 
of  the  steady  state  inflow  rate,  Q,  and  the  saturated  hydraulic  conductivity,  Ks,  yielded  a 
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linear  relationship.  Comparison  of  inverse  solution  results  for  Kg  from  the  laboratory  aquifer 
and  field  tests  of  the  cone  permeameter  with  Ks  values  predicted  with  this  relationship 
differed  by  less  than  a  factor  of  three.  Sensitivity  analyses  showed  that  this  relationship  was 
significantly  influenced  by  the  van  Genuchten  parameters  a  and  n  while,  but  not  by  the 
values  of  Of  and  ds-  (Al-Houri,  1999).  More  work  is  needed  on  this  topic  to  determine  if  a 
reasonably  accurate  empirical  analysis  method  can  be  developed.  A  rapid  assessment 
technique  would  increase  the  general  usefulness  of  the  cone  permeameter  by  decreasing  the 
computational  effort  required  to  analyze  the  data  when  only  knowledge  of  Kg  is  needed. 

•  The  results  of  the  project  show  that  the  cone  permeameter  method  can  be  used  to  predict  the 
soil  hydraulic  properties  in  sandy  soils.  We  have  begun  field  testing  the  prototype  and 
analysis  method  in  a  loamy  soil  to  test  the  method  in  soils  with  greater  fines  contents. 
Confirmatory  laboratory  and  field  tests  are  ongoing. 


Figure  3.  Cone  permeameter  set  up  for  field  testing. 
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Figure  4.  Observed  and  simulated  responses  for  cone  permeameter  Test  B  at  Site  1  in  Poinsett 
State  Park  (Kodesova  et  al.,  1999b). 
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Figure  5.  Soil-moisture  retention  curves  resulting  from  standard  methods  and  inversion  of  flow 
data  from  five  cone  permeameter  tests  at  Site  1  in  Poinsett  State  Park.  Curves  A,  C,  and  D  are 
from  analysis  of  the  infiltration  part  of  the  experiments,  while  B  and  E  are  from  both  infiltration 
and  redistribution  (dB  and  dE  are  the  drying  branches).  Curves  CR  and  MO  are  from  capillary 
rise  and  multi-step  outflow  experiments,  respectively  (Kodesova  et  al.,  1999b). 
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Table  1.  Saturated  hydraulic  conductivity  values  obtained  via  inversion  of  cone  permeameter 
data  and  standard  test  methods  (ho  =  supply  pressure  head  value). 


Test  Method 

Ks  (cm/sec) 

Multi-Step  Outflow  Tests  (2  samples) 

0.0012-0.0027 

Guelph  Permeameter  Tests  (6  locations) 

0.0025-0.0039 

Laboratory  Falling  Head  Tests  (9  samples) 

0.0013-0.0044 

Cone  Permeameter  Test  A:  ho=  30,  50  cm 

0.0016 

Cone  Permeameter  Test  B:  ho=  30,  50  cm,  &  redistribution 

0.0016 

Cone  Permeameter  Test  C:  ho=  30,  50  cm 

0.0011 

Cone  Permeameter  Test  D:  ho=  21, 108  cm 

0.0036 

Cone  Permeameter  Test  E:  ho  =  21,  80  cm,  &  redistribution 

0.0010 

Army  Relevance 

The  products  of  this  research  project  have  the  potential  to  enhance  the  Army's  ability  to 
characterize  and  remediate  contaminated  sites.  When  this  technology  fully  developed,  it  will 
have  numerous  advantages  over  existing  methods  for  soil  hydraulic  property  measurement.  As 
with  other  CPT  tools,  the  cone  permeameter  is  minimally  intrusive  and  does  not  result  in  the 
removal  of  potentially  contaminated  materials.  This  feature  safeguards  the  welfare  of  operating 
crews  and  reduces  site  investigation  costs.  It  is  applicable  to  depths  of  30  m  or  more  below 
ground  surface,  depending  on  soil  resistance  to  penetration.  Since  the  analysis  is  based  on 
transient  flow  data,  time  and  water  requirements  are  minimized.  Finally,  the  inverse  solution 
approach  makes  it  possible  to  simultaneously  determine  the  soil-moisture  retention  and  hydraulic 
conductivity  curves  in  unsaturated  soils,  something  not  currently  possible  with  other  methods. 
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